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I. INTRODUCTION 



The phenomenology of QCD exclusive processes depends in a crucial way on the analytic 
properties of hadronic (hard) scattering amplitudes as functions of the strong running cou- 
pling. A perturbatively calculable short-distance part of the reaction amplitude at the parton 
level is isolated either by subtraction or by factorization. To get a quantitative interpretation 
of such quantities in practice and compare them with experimental data, one has to get rid 
of the artificial Landau singularity at Q 2 = Aq CD (Aqcd = A in the following), where Q 2 is 
the large mass scale in the process. A proposal to solve this problem (in the spacelike region) 
without introducing exogenous infrared (IR) regulators, like an effective, or a dynamically 
generated, gluon mass [1] (see, for instance, [2, 3, 4, 5, 6, 7, 8, 9, 10] for such applications), 
was made by Shirkov and Solovtsov (SS) [11, 12, 13], based on general principles of local 
Quantum Field Theory. This theoretical framework — termed Analytic Perturbation Theory 
(APT) — was further expanded beyond the one-loop level of two-point functions to define 
an analytic 1 coupling and its powers in the timelike region [14, 15, 16, 17, 18, 19, 20, 21], 
embracing previous attempts [22, 23, 24, 25, 26, 27] in this direction. 2 

However, first applications [30, 31] of this sort of approach to three-point functions, 
beyond the leading order of QCD perturbation theory, have made it clear that, ultimately, 
there must be an extension of this formalism from the level of the running coupling and 
its powers to the level of amplitudes. The reason is that in three-point functions at the 
next-to-leading order (NLO) level, and beyond, logarithms of a distinct scale (serving as 
the factorization or evolution scale) appear that though they do not change the nature of 
the Landau pole, they affect the discontinuity across the cut along the negative real axis 
— oo < Q 2 < 0. On account of factorization, we expect that this effect should be small, of 
the order of a few percent, because any change caused by the variation of the factorization 
scale should be of the next higher order. However, to achieve a high-precision theoretical 
prediction, one should reduce this uncertainty, lifting the limitations imposed by the lack 
of knowledge about uncalculated higher-order corrections. To encompass such logarithmic 
terms in the "analytization" procedure, one should demand the analyticity of the partonic 
amplitude as a whole [32, 33] and calculate the dispersive image of the coupling (or of its 
powers) in conjunction with these logarithms. This Karanikas-Stefanis (KS) "analytization" 
scheme effectively amounts to the generalization of APT to non-integer powers of the running 
coupling: Fractional APT (FAPT), as we shall show below. 

In this work we expand the Shirkov-Solovtsov "analytization" approach to include the 
dispersive images of such terms, using as a case study the pion form factor at NLO in the MS 
scheme with various renormalization-scale settings and also in the ay-scheme [34]. To this 
end, we contrast the KS "analytization" with the naive [30, 31] and the maximal [35] "analy- 
tization" procedures and work out their key mutual differences as they first appear in NLO, 
while a fully-fledged analysis of FAPT is given in an accompanying paper [36]. We argue 
that augmenting the MS scheme with the KS "analytization" prescription provides an opti- 
mized method to calculate perturbatively higher-order corrections to partonic "observables" 
in QCD because it practically eliminates all scheme and scale-setting ambiguities owing to 
the renormalization and factorization scales. It is worth emphasizing at this point that the 
focus of Ref. [32] was on the calculation of power corrections to the pion's electromagnetic 



1 The term 'analyticity' is used here as a synonym for 'spectrality' and 'causality' [12]. 

2 A somewhat different approach was reviewed recently in [28] ; see also [29] . 
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form factor. Such contributions are outside the scope of the present investigation. 

The plan of this paper is as follows. In Sec. II we review the convolution formalism for the 
calculation of the short-distance part of the pion form factor within perturbative QCD at 
NLO. In Sec. Ill we discuss the Shirkov-Solovtsov type "analytization" procedures [11, 30, 
31, 32, 33, 35] and work out their mutual differences, focusing on the KS "analytization" and 
its properties. This discussion extends and generalizes the original KS analysis that covered 
only the LO of the perturbative expansion of the pion form factor and ignoring evolution. 
Section IV contains the results for the factorized pion form factor in different schemes and 
with different scale settings, employing the KS "analytization" in comparison with those 
based on APT and also standard QCD perturbation theory in NLO. Our conclusions with 
a summary of our main results are presented in Sec. V. Some important technical details 
are collected in three appendices. 



II. FACTORIZABLE PART OF THE PION FORM FACTOR AT NLO IN STAN- 
DARD QCD PERTURBATION THEORY 

The leading-twist factorizable part of the electromagnetic pion form factor can be ex- 
pressed as a convolution in the form [37, 38] 

if act (Q 2 ; /4) = i4) ® T H (x, y, Q 2 ; & /4) ® * K (y, , (2.1) 

x y 

where ® denotes the usual convolution symbol (A(z)®B(z) = dzA{z)B{z)) over the 

longitudinal momentum fraction variable x (y) and /i F represents the factorization scale at 
which the separation between the long- (small transverse momentum) and short-distance 
(large transverse momentum) dynamics takes place, with /xr labelling the renormalization 
(coupling constant) scale. The nonperturbative input is encoded in the pion distribution 
amplitude (DA) $ n (y, whereas the short- distance interactions are represented by the 
hard-scattering amplitude T H (x, y, Q 2 ; /i R ). This is the amplitude for a collinear valence 
quark- ant iquark pair with total momentum P struck by a virtual photon with momentum q, 
satisfying q 2 = —Q 2 , to end up again in a configuration of a parallel valence quark-antiquark 
pair with momentum P' = P + q. It can be calculated perturbatively in the form of a power- 
series expansion in the QCD coupling, the latter to be evaluated at the reference scale of 
renormalization // R : 

T H NLO (x, y, Q 2 - & /4) = a s (&) T<°> (x, y, Q 2 ) + r a> (a . y> Q 2. ^ ^ (2 2) 

The leading-order (LO) contribution to Tu(x,y,Q 2 ; /ip) reads 

Tl°\ X ,y,Q-l = ^± = ±tf( X ,y), (2.3) 



where 
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C F = (iV c 2 - 1) /2N C = 4/3, C A = N c = 3 are the color factors of SU(3) C , and the notation 
z = 1 — z has been used. The usual color decomposition of the NLO correction [39] — marked 
by self-explainable labels — is given by (omitting the variables x and y) 



QX ] (Q 2 ; A, 4) = c F ig" 4 1 ' 91 (^) + c 4''°' , 



(2.5) 



where Cq = (Cf — Ca/2) and &o is the first coefficient of the /3 function, see Appendix 
A, Eq. (Al). Here we explicitly factorized out a trivial 1/Q 2 dependence and used for the 
coefficients in front of each factor the notation t H with appropriate superscripts. 

With reference to the application of the Brodsky-Lepage-Mackenzie (BLM) [40] scale 
setting in fixing the renormalization point later on, we single out the &o-P ro P°rtional (i.e., 
the iV/-dependent) term, given by 



with 



(0) 



--\n(xy) 



and present the color singlet part of tu in the form 



*(1,F) 
y H,2 



/4 
Q 2 



4 0) (*,!/) 



2(Wln(xy)) In 



/4 



(2.6a) 

(2.6b) 
(2.6c) 

(2.7a) 
(2.7b) 



Explicit expressions for t^'^ (x,y) and for the color non-singlet part, t^' G \x,y), cf. Eq. 
(2.5), are supplied in Appendix B (see Eqs. (Bl), (B2)). 

The scaled hard-scattering amplitude, Eq. (2.2), truncated at the NLO and evaluated at 
the renormalization scale /x^ = ArQ 2 , reads 

Q 2 T H NLO (x, y, Q V 2 , XrQ 2 ) = a s (A R Q 2 ) i 0) (x, y) + C F f £f (x, y; ^ 



+ 



a 2 s (XrQ 2 
An 



{b t^\x,y;\ R )+t ( £ G \x,y)} , (2.8) 



where we have introduced the shorthand notation 

t (FG)/ \ _ ^ ,(1,F), 



f^'Or, y) = C F #f (x, y) + C G ^(x, y) . 



(2.9) 



To calculate the factorizable part of the pion form factor, one has to convolute this expression 
with the pion DA for each hadron in the initial and final state. In leading twist 2, the pion 
DA at the normalization scale ^ ~ 1 GeV 2 is given by 

y^(x,/i 2 ) = 6x(l - x) [l + a 2 (/i 2 ) C 2 3/2 (2x - 1) + a 4 (/i 2 ) C 4 3/2 (2:r - 1) + . . .] , (2.10) 



with all nonperturbative information being encapsulated in the Gegenbauer coefficients a n . 
In this analysis we use those coefficients determined before by Bakulev, Mikhailov, and 
Stefanis (BMS) in [41] with the aid of QCD sum rules with nonlocal condensates: 

a* MS = 0.20, a? MS = -0.14, a™ s = 0,n>4, (2.11) 

where the vacuum quark virtuality A 2 = 0.4 GeV 2 has been used. This set of values was 
found [42, 43] to be consistent at the la level with the high-precision CLEO data [44] on the 
pion-photon transition form factor, with all other model DAs being outside — at least — the 
2a error ellipse (see [45] for the latest compilation of models in comparison with the CLEO 
and CELLO [46] data). Notice that the particular parameterization (shape) of the pion DA 
chosen is irrelevant for the considerations to follow. 



III. ANALYTICITY OF PARTONIC AMPLITUDES BEYOND LO 
A. Analytic Running Coupling in QCD 

The main stumbling block in applying fixed-order perturbation theory at low momenta 
Q 2 is the non-physical Landau singularity of the running strong coupling at Q 2 = A 2 , which 
entails the appearance of IR renormalons in the perturbative expansion. To ensure the 
analyticity of the coupling in the infrared, one can follow different strategies [11, 25, 27, 
47, 48, 49, 50, 51] all based on the basic assumption that the physical coupling should stay 
IR finite and analytic in the whole momentum range, though its precise value at Q 2 = is 
still a matter of debate [11, 21, 28, 52, 53]. Imposing the analyticity of the coupling in the 
sense of Shirkov and Solovtsov [11], we replace the strong running coupling and its powers 
by their analytic versions: 



i 



Im [/(-<r)] 



da , 



n Jo a + Q 2 -ie 

where the loop order is explicitly indicated by the superscript n in parenthesis and 

1 A 2 



4" (<2 2 ) = £ 

°0 



+ 



ln(Q 2 /A 2 ) A 2 -Q 2 



a. 



(Q 2 ), 



(3.1) 



(3.2) 



with the last step connecting to the SS notation [11], and a s (0) = 47r/&o- The two-loop 
running coupling in standard QCD perturbation theory can be expressed [19] in terms of 
the Lambert function W_x to read 



a 



(2) 



(Q 2 ) 



47T 



(3.3) 



For some more explanations we refer the interested reader to [42], Appendix C, Eqs. (C15) 
and (C20). Then, the analytic image of the fcth power of the coupling [14] is obtained from 
the dispersion relation 



pIV) 

+ Q 2 - ie 



(3.4) 
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TABLE I: Parameters entering Eqs. (3.6) and (3.7) for the value Aq£ d 3 = 400 MeV. 



Parameters 


Jit 
c 21 


A 2 i 


Jit 
L 22 


A 22 


Values 


-1.015 


67 MeV 


-1.544 


34.5 MeV 



with the spectral density 



k / , \ k 



In the numerical calculations below, we use an approximate form suggested in [35]: 

,(2,fit) (Q2) = 4vr f 1 1 1 

1 lW j &o 1 ^ [In (Q 2 / A ii) , 4\] + 1 - exp(£ [In (Q 2 / A 2i) , 4\] ) J ' ^ 



A VM(QZ) = f 47r V( 1 exp(£[ln(Q 2 /Al 2 ),c^]) 

2 V / U fin (Q 2 /AL) , c<*l 2 [l - exp(£ fin (Q 2 /A 2 ) , c*l 



(3.7) 



[In (Q 2 /A 2 2 ) , c*] [1 - exp(£ [in (Q 2 /A 2 2 ) , <*] )] ' 
where the values of the fit parameters are listed in Table I and 

£[L,c] = L + c\nVL 2 + Att 2 . (3.8) 



B. "Analytization" Procedures 

Let us now see how analyticity can be implemented on the parton-level pion form factor 
in NLO accuracy of perturbative QCD. We discuss three "analytization" procedures: 3 

• Naive "analytization" [30, 31, 35] 

[Q 2 T H (s,y,g 2 ;/4, A R Q 2 )]~ = A?\\ R Q 2 ) t£\x, y) 



+ 



4:71 



b t^\x, y;\ R )+ 4 FG) (x, y) + C F fgf (x, y; ^ 



(3.9) 



• Maximal "analytization" [35] 



4 2) (A R g 2 

4tt 



ftb t ( i' P) (x.yy:A H ) + /[r T '(.«--/y) + CV^o' (•'-.//: Q2 



f (FG), 



,(1,F) 



.(3.10) 



3 One should not worry about the factor l/Q 2 because under "analytization" it reproduces itself, i.e., 

[[/(0 2 )] an /Q 2 ] an =[/(Q 2 )] an /0 2 - 
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• Amplitude "analytization" proposed by Karanikas and Stefanis in [32, 33]. 

The first method replaces a s and its powers by the Shirkov-Solovtsov analytic coupling [1 1] 
and its powers, whereas the second one uses for the powers of a s their own analytic images, 
transforming this way the power-series expansion in [a s {Q 2 )] n in a functional expansion 
in terms of the functions A n (Q 2 ) [13, 16]. Imposing analyticity in the sense of Karanikas- 
Stefanis [32] , differs from the previous two approaches in that it demands the whole partonic 
amplitude has the correct analytical behavior as a function of Q 2 . This entails the "analyti- 
zation" of terms of the form (Q 2 )j ln(Q 2 //i F ), which appear in exclusive amplitudes 

at NLO of QCD perturbation theory and contain an additional scale, //p. There are, in 
principle, two possibilities how to proceed any further. One option is provided by setting 
/xp ~ Q 2 and then face the problem of "analytization" of terms like [a s (Q 2 ) /ct s (/x 2 l )]' ? , where 
f] = 7^/(260) is a fractional number, as discussed in [36]. Another possibility is to fix the 
factorization scale /xp at some value and then to redefine the original Shirkov-Solovtsov "an- 
alytization" procedure in order to take the dispersive image of the coupling (or of its powers) 
together with these logarithmic terms. This second route is followed in the present work. 
It is important to note that the KS "analytization" procedure reduces in LO of fixed-order 
perturbation theory to the maximal one, as shown in [32], provided evolution effects of the 
pion distribution amplitudes are ignored. 

Applying now this generalized "analytization" concept, we get 



[Q 2 T H (x, y, Q 2 ; & XrQ 2 )}™ = A?\\ R Q 2 ) t^(x, y) 



+ 



4 2) (ArQ 2 ) 



47T 



(b t^\x,y;\ R )+t ( i G \x,y)) 



+ 



(2) 



(ArQ 2 ))' 



4TX 



C F t$\x,y) (6 + 21n(xy))ln^ 



(3.11) 



KS 



In order to have the same scale argument in the logarithmic term as in the running coupling, 
we substitute ln(Q 2 //ip) = \n(\ R Q 2 /A 2 ) — ln(A R //p / A 2 ) to obtain 



[Q 2 T H (x, y, Q 2 ; /4, A R Q 2 )] £ = A? (A R Q 2 ) C v) 



f(°)/ 



+ 



+ 



4 2) (a r q 2 ) 



4tt 



b 4 1>/J) (x, y; A R ) + t£ ^> (x, y) - C F t%> (x, y)(6 + 2 \n(xy)) In 



f (FG) 



(0), 



A R /ip 



(2) 
Oil 



(A R g 2 



(°), 



4lX 



C F t H 



x, y) (6 + 2 ln(x y )) In 



ArQ 2 

A 2 



J KS 



Finally, we arrive at 
[Q 2 T H (x,y,Q 2 ;nl, \ R Q' 



an 
KS 



+ 



+ 



4 2) (a r q 2 



4lX 



A?\\ R Q 2 )4\x,y) 

b (x, y; Ar) + 4 FG) (x, y) + C F fgf 



Ag (A R Q 2 
4tt 



C F t£\x,y) (6 + 21n(xj/)) 



(3.12) 



/4 
Q 2 



(3.13) 



7 



0.5 
0.4 
0.3 

0.2 
0.1 

■ 



4 2) (Q 2 



0.5 
. 4 
. 3 
0.2 
. 1 




Af (Q 2 



1 - (2) , n 2^ 
\ 

\ 

\ 

• \ 



(b) 



10 



20 



30 



4 



[GeV 2 



50 



10 



20 



30 



40 



50 



Q 2 [GeV 2 



FIG. 1: (a): Results for the analyticized logarithmic term, C 2 (Q 2 )' usm g three different "analytization" 



procedures: the one-loop approximate KS "analytization", C 2 



(1) approx 



loop approximate KS "analytization" , £ 2 ' 



(2) approx 



(Q 2 ) (red dash-dotted line), the two- 



(Q 2 ) (green dashed line), and the exact two-loop BMKS 



"analytization" (Q 2 ) (blue solid line). For comparison, we also show here the corresponding maximal 
"analytization" curve (dotted line), (b): Results are shown for the corresponding analyticized contributions 
A 2 2 ^ {Q 2 )> using the same three "analytization" procedures for (Q 2 ) a s in panel (a). 



where the deviation from Eq. (3.10) is encoded in the term 



A< 2 > (Q 2 ) = Cf (Q 2 ) - A? (Q 2 ) In [Q 2 /A 2 ] 



with 



6 2) 



(Q 2 ) 



W 2) [Q 2 )Y^[% 



J KS 



An 



a 



(2) 



(Q 2 



21 an 



a 



(i) 



(Q 2 ) 



(3.14) 



(3.15) 



J KS 



It is important to distinguish between the two contributions in Eq. (3.14). The first amounts 
to the "analytization" of the product of the coupling with a logarithm, or equivalently 
of fractional powers of the coupling, as shown in [36]. The second bears an additional 
logarithmic dependence on the momentum scale Q 2 relative to the expression obtained with 
the maximal "analytization" procedure. The subscript KS in the last equation signifies 
that this expression should be analyticized according to the KS prescription. To obtain a 
clearer idea of its meaning and demonstrate its essence, the "analytization" is performed in 
three incremental steps. First, a simplified version of this expression is considered, which 
results by provisionally replacing the two-loop coupling in the numerator by its one-loop 
counterpart. Then, the ratio of the couplings after "analytization" reduces to (dash-dotted 
line in Fig. la) 



p (1) approx 



An 



(3.16) 



Second, we discuss an analogous situation, in which the one-loop coupling in the denom- 
inator is (inconsistently) traded for its two-loop counterpart. In this case, the ratio of the 
couplings after "analytization" becomes (dashed line in Fig. la) 



(2) approx 



Air 



(Q 2 ) = yA 2 > (Q 2 ). 



(3.17) 



S 



Finally, we provide the exact result for the KS "analytization" of expression (3.15) (solid 
line in Fig. la), with the derivation presented in Appendix C, while more general expressions 
are given in [36]: 



4 2) (Q 2 ) 



An 



A 2) (Q 2 ) +cij^fc (Q 2 ) 



(3.18) 



where 



u m = e 



n>0 



^(2)C(-n-l)- 



d((-n - 1) 
dn 



[~ln(g 2 /A 2 )] r 

r(n + i) 



(3.19) 



and ((z) is the Riemann zeta-function. Equation (3.14) is illustrated in Fig. lb for the 
different expressions of £ 2 {Q 2 ) given by Eqs. (3.16), (3.17), and (3.18), using the same line 
designations as in Fig. la. Let us close this discussion by commenting that in the region 
where there are experimental data available [54, 55] (i.e., well below 10 GeV 2 ), Eq. (3.14) is 
governed by {Q 2 ), which entails a small enhancement of the hard-scattering amplitude 
for Q 2 < 7.25 GeV 2 . 



IV. FACTORIZED PION FORM FACTOR AT NLO STANDARD AND ANA- 
LYTICIZED 

The calculation of the factorized pion form factor proceeds in terms of Eq. (2.1) and 
involves the convolution of expression (3.10) for the maximal "analytization" case, or ex- 
pression (3.13) for the KS "analytization" case with the pion DA for which we employ in 
both cases the BMS parameterization [41], as discussed in Sec. II. On that basis, we can 
obtain the scaled, factorized part of the pion form factor, Q 2 F^ act (Q 2 ;^ = XrQ 2 ), using 
Eq. (2.8) and the following set of substitutions: 4 

8tt/ 2 (l + a 2 + a 4 ) 2 ; (4.1) 

8 7T f 2 (1 + a 2 + a 4 ) [3 + (43/6)a 2 + (136/15)a 4 ] ; (4.2) 
8 7T f 2 [ - 15.67 - a 2 (21.52 - 6.22 a 2 ) 

- a 4 (7.37 - 37.40 a 2 - 33.61 a 4 ) ] . (4.3) 

Notice that evolving the BMS pion DA from the initial scale /i 2 , to the scale /j| at the 
NLO level will generate higher Gegenbauer harmonics of the form xxC^Qx — 1) with 
n > 3. However, we have shown in [35] (see also [56]) that for the calculation of the pion 
form factor it is actually sufficient to restrict ourselves to the LO evolution and neglect NLO 
evolution effects. Hence, for our purposes in the present analysis, we set 

7i°V(2M 

(4.4) 



4 Here, we write for the sake of brevity a 2 = 02 MS (^f) an d a 4 — a 4 MS (Mp) an d use the values given in Eq. 
(2.11). 



4 FG, (^) 



fl2n(/4) = «2n(/io) 
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The lowest-order anomalous dimensions can be represented in closed form by 

2 



7i 0) = 2C, 



451(71 + 1) -3 



(n + l)(n + 2) 



(4.5) 



with Si(n + 1) = Y^=i V* = + 2) — VK-Qj while the function is defined as ip(z) = 
dlnT{z)/dz. 

Following the master plan for "analytization" , exposed in the previous section, we obtain 
the following expressions for the factorized pion form factor: 

• Naive "analytization" [30, 31, 35]: 

K act (Q 2 ;A R g 2 )] NaivAn = 4 2) (a r q 2 )^ lo (q 2 ) 



+ 



7T 



4 2) (A R Q 2 ) ~^ LO (Q 2 ,/4;A R ) . (4.6) 



Maximal "analytization" [35]: 

K act (Q 2 ;A R Q 2 )] MaxAn = 4 2) (ArQ 2 )^ lo (Q 2 ) 



+ -4 2) (A R Q 2 )^ NLO (Q 2 ,/4;A R ) 



7T 



(4.7) 



KS amplitude "analytization" (this work) — cf. Eqs. (3.13) and (3.15): 

K act (Q 2 ;A R Q 2 )] KS = 4 2) (ArQ 2 )^ lo (Q 2 ) 



+ -4 2) (ArQ 2 )^ nlo (Q 2 ,/4;A r ) 

71 

+ ( X *Q 2 ) Af jtNLO /q2) 

7T 



Here we use the following notations: 



^° (Q 2 ) = ^ (1 + a 2 + a 4 ) 2 ; 



(4.8) 



(4.9) 



^ NLO (Q 2 ,/i|;A R ) 



2tt/ 2 

Q 2 



6 (1 + a 2 + a 4 ) 2 (In A R - In A B lm(o2, a 4 )) - 15.67 



- a 2 (21.52 - 6.22 a 2 ) - a 4 (7.37 - 37.40 a 2 - 33.61 a 4 ) 

2 



+ A F ^ NLO (g 2 )ln% 



(4.10) 



and we explicitly display the contribution due to t^'p (x,y; Hp/Q 2 ), see Eq. (2.7b): 



A F ^ LO (Q 2 ) = ~^Sc F (l + a 2 + a 4 )[(25/3)a 2 + (182/15)a 4 ] . (1.11) 



Q 2 

In order to make our formulas more compact, we implement the BLM scale: 

5 3 + (43/6)a 2 + (136/15)a 4 " 



ABLM(a2,a 4 ) = exp 



1 + a 2 + a 4 



(4.12) 
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FIG. 2: Results for the ratio of the factorized pion form factors, using two different "analytization" proce- 
dures: KS "analytization" and "maximal analytization" , F^ s (Q 2 ) /.pMax-an (Q 2 ). The designations are: red 
dash-dotted line — one-loop approximation of the KS logarithmic term (C^ a PP rox (Q 2 )); green dashed line — 
two-loop approximation (£( 2 ) a PP lox (Q 2 )); blue solid line — exact two-loop KS "analytization" (£2 (Q 2 ))- 
Left panel: default scale setting (A = 1); middle panel: BLM scale setting; right panel: ay scheme. The 
factorization scale fi F is set equal to 5.76 GeV 2 [57]. 



The "analytization" augmented perturbation theory works very well. This is illustrated 
by the results in Figs. 2, 3, and 4. The first of these figures compares the specific issues of 
the KS "analytization" procedure relative to those of the maximal one for the ratio of the 
corresponding factorized form factors. A few words are in order here. One sees that using 
the default MS scheme, the KS "analytization" procedure yields a result almost coincident 
with that provided by the maximal one. On the other hand, in the BLM scheme and also 
in the ay scheme, the KS prediction is smaller by a few percent. Moreover, one observes 
by comparison with Fig. 11, right panel in Ref. [35] that the BLM prediction, which in 
the maximal procedure was the largest one, becomes in the case of the KS prescription 
comparable with the prediction of the default scheme. As a result, the inherent theoretical 
uncertainties due to the involved perturbative parameters, defining a renormalization scheme 
and scale setting, are further reduced. A second important feature of the KS procedure is 
that the dependence of F^ act (Q 2 ) on the factorization scale is almost diminished, as indicated 
in Fig. 3. Indeed, varying the factorization scale from 1 GeV 2 to 10 GeV 2 , the form factor 
changes by a mere 1.5 percent. Even setting the factorization scale to the theoretical value 




FIG. 3: Results for the factorized pion form factors, using two different "analytization" procedures: KS 
"analytization" (blue solid line) and "maximal analytization" (dotted red line) for different values of the 
factorization scale: in the left panel fi F = 1 GeV 2 , in the middle one-/Xp = 1 GeV 2 , and in the right 
one-/z| = 10 GeV 2 . For all panels we show results for the BLM scale setting. 
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of 50 GeV 2 , the induced variation in the form-factor magnitude reaches just the level of 
about 2.5 percent. In the case of the maximal "analytization" procedure, the dependence 
on the factorization scale is also a mild one, but the corresponding variation is, in round 
terms, two times larger. 

The fourth figure demonstrates the impact of "analytization" on the factorized pion's 
electromagnetic form factor, using various "analytization" prescriptions. The dashed line 
denotes the prediction obtained with standard QCD perturbation theory in the MS scheme 
and applying the default scale setting n\ = Q 2 . The naive "analytization" prediction is 
represented by the dash-dotted line and the analogous one for the maximal "analytization" 
by the solid line below it. The result of the calculation according to the KS "analytization" 
practically coincides with that of the maximal one. This behavior is also reflected in Fig. 
2, where we see that the differences among the three "analytization" procedures are of the 
order of a few percent in the whole Q 2 range considered. 

Note that as regards the whole pion form factor, i.e., taking into account also the soft 
part, the differences would be further reduced. For full details the reader is referred to [35]. 



. 35 
0.3 

0.25 
. 2 



2 3 4 5 6 7 8 

Q 2 [GeV 2 ] 

FIG. 4: Results for the factorized pion form factor, scaled with Q 2 , and assuming the default scale setting 
(Mr = Q 2 ) m standard perturbation theory and APT. The latter is implemented in terms of two different 
"analytization" procedures: naive "analytization" and maximal "analytization". The designations are: blue 
dashed line — standard perturbation theory; green dash-dotted line — naive APT; red solid line — maximal 
APT. The prediction obtained with the KS "analytization" is too close to that found with the maximal one 
to differentiate these curves graphically. The factorization scale /i| is set equal to 5.76 GeV 2 . 



V. SUMMARY AND CONCLUSIONS 

We have discussed different "analytization" procedures to ensure the analyticity of the 
factorized electromagnetic pion form factor at NLO of QCD perturbation theory. The main 
features and relative merits of each "analytization" concept following from the presented 
analysis are: 

• The naive "analytization [30, 31] retains the power-series expansion of perturbative 
QCD, but replaces (ai n) (Q 2 )) m by (A^ (Q 2 )) m . As it was shown in [30, 31], this 
reduces the value of the NLO correction, though the sensitivity to the renormaliza- 
tion scheme adopted and the renormalization scale-setting chosen is still substantial, 
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resulting into a rather strong variation of the form-factor predictions [35]. Moreover, 
this procedure does not respect nonlinear relations of the coupling because these cor- 
respond to different dispersive images. 

• The maximal "analytization" [35] trades the power-series expansion for a functional 
non-power-series expansion m terms of A { m\Q 2 ) [11, 16, 17], minimizing the variation 
of the form-factor predictions owing to the renormalization scheme and scale setting. 
It is, however, insufficient to cure logarithms of the momentum scale multiplying the 
running coupling. Such terms modify the spectral density, i.e., the discontinuity across 
the cut along the negative real axis and have therefore to be taken into account. 

• Applying the "analytization" procedure at the level of the partonic amplitude itself 
[32, 33], bears all advantages of the maximal "analytization" plus a further reduced 
dependence on the perturbative scales — especially the dependence on the factoriza- 
tion scale. This has been verified by explicit calculation. We have employed the MS 
scheme with various scale settings and also the ay scheme. In addition, we have var- 
ied the factorization scale in the range 1 — 10 GeV 2 . While the predictions for the 
factorized pion form factor, calculated with the maximal procedure, were affected by 
this variation on the level of 3%, their counterparts, derived with the KS prescription, 
were influenced by less than 1%. Though the KS method does not really "gain up" 
relative to the maximal "analytization" procedure with respect to the factorized pion 
form factor, as one observes from Fig. 4, it is able to further improve the perturba- 
tive treatment because it extends the notion of analyticity to non-integer powers of 
the strong running coupling — FAPT. Such powers become relevant when one has to 
calculate the analytic image of powers of the strong coupling in combination with 
logarithms, the latter first appearing at NLO of fixed-order perturbation theory, or in 
terms of evolution factors [36]. Hence, the KS "analytization" requirement treats all 
logarithms that have a non-zero spectral density, and hence modify the discontinuity 
across the cut along the negative real axis, on the same footing and irrespective of 
their source being it the running coupling (and its powers), or logarithms entailed by 
ERBL or DGLAP evolution. 

In conclusion, the KS "analytization" enables the variation of the factorization scale and 
the choice of various renormalization schemes and scale settings, including the BLM one, with 
undiminished quality of the theoretical predictions from scheme (scale) to scheme (scale), 
virtually eliminating the dependence on such parameters and upgrading the MS scheme to 
an optimized factorization and renormalization scheme. From a broader perspective one 
may interpret these findings as indicating that the analyticity of the partonic three-point 
function is as important and fundamental as the underlying symmetries of the theory and 
should be preserved together with them in the maximal possible way. 
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APPENDIX A: QCD (5 FUNCTION AT NLO 



The first coefficients of the (3 function are 



11 4 

&o = yC A -3T R iV / , 



34 



AC F + ^-C A \T R N f . 



(Al) 



Here, T R = 1/2 and Nf denotes the number of flavors, whereas the expansion of the (3- 
function in the NLO approximation is given by 



47T 



47T 



(A2) 



APPENDIX B: NLO CORRECTION TO THE PION FORM FACTOR 



Here we present the detailed expressions for the color decomposition of the NLO correction 
to the hard amplitude T H , which describe the factorized part of the pion form factor [35, 39] 
(see Eqs. (2.5)-(2.7b)): 
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The functions H(x,y) and R(x,y) are defined by 
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and 
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APPENDIX C: "ANALYTIZATION" OF POWERS OF THE COUPLING MUL- 
TIPLIED BY LOGARITHMS 

We present here the derivation of {Q 2 ), done in collaboration with S. Mikhailov. To 
this end, let us first introduce 

a s (Q 2 )=^a s (Q 2 ) . (CI) 
14 



For this quantity we can write a renormalization group solution in the form 



A 2 



1 + Cl ai 2) (Q 2 ) 



(C2) 



Expanding the expression ln[l + Cia^ (Q 2 )] and retaining terms up to order a 2 s , we find 



K 2) (Q 2 )} 2 ^ (jp) = ^ (Q 2 ) + K 2) (Q 2 )} 2 <> 1- K 2) (Q 2 )} 

To get rid of the logarithm, we use the following trick 



(C3) 



K 2) (Q 2 )T^(%)= a( s ) m + ^l'P (Q 2 )] 



2+e 



and return to the original coupling to obtain 



de 



e=0 



(C4) 



Now we can proceed with the "analytization " of the term [af^ (<3 2 )] 2 ln (Q 2 ), giving rise to 
analytic expressions for non-integer powers of the coupling, i.e., 



{[^(^] a h(g)}=£>iP> (<?)+«. 



(C6) 



£=0 



Using the representation [36] 



[-ln(Q 2 /A 2 )]" 
T(n + 1) 



n>0 



and performing the differentiation, we finally obtain 



with / £ (Q 2 ) being defined in Eq. (3.19). 
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